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Accretion and diffusion in white dwarfs 

New diffusion timescales and applications to GD362 and G 29-38 

D. Koester 

Institut fiir Theoretische Physik und Astrophysik, University of Kiel, D-24098 Kiel, Germany 

ABSTRACT 

Context. A number of cool white dwarfs with metal traces, of spectral types DAZ, DBZ, and DZ have been found to exhibit infrared 
excess radiation due to circumstellar dust. The origin of this dust is possibly a tidally disrupted asteroid that formed a debris disk now 
supplying the matter accreting onto the white dwarf. To reach any clear conclusions from the observed composition of the white dwarf 
atmosphere to that of the circumstellar matter, we need a detailed understanding of the accretion and diffusion process, in particular 
the diffusion timescales. 

Aims. We aim to provide data for a wide range of white dwarf parameters and all possible observed chemical elements. 
Methods. Starting from atmosphere models, we calculate the structure of the outer envelopes, obtaining the depth of the convection 
zone and the physical parameters at the lower boundary. These parameters are used to calculate the diffusion velocities using calcula- 
tions of diffusion coefficients available in the literature. 

Results. With a simple example, we demonstrate that the observed element abundances are not identical to the accreted abundances. 
Reliable conclusions are possible only if we know or can assume that the star has reached a steady state between accretion and diffu- 
sion. In this case, most element abundances differ only by factors in the range 2-4 between atmospheric values and the circumstellar 
matter. Knowing the diffusion timescales, we can also accurately relate the accreted abundances to the observed ones. If accretion has 
stopped, or if the rates vary by large amounts, we cannot determine the composition of the accreted matter with any certainty. 

Key words, stars: white dwarfs - stars: abundances - accretion - diffusion 



1. Introduction 

Heavy elements in cool white dwarfs 25 000 K) should dif- 
fuse downward in the atmospheres due to gravitational settling 
in the high gravitational fields (Schatzman, 1945). Nevertheless, 
one of the first three "classical" white dwarfs, van Maanen 2, 
has been found to exhibit strong Can resonance lines. This 
star has a helium-rich atmosphere and is the prototype of the 
spectral class DZ. Although an example of a star containing a 
hydrogen-rich atmosphere with metal traces, G74-7, had been 
known for a consid erable period of time (see the analysis by 
iBilleres et all 1 19971) . it was not until 19 97, that DAZ \ yas fi- 
nally established as an important class jKoester et all 119971; 
iHolberget all 119971) . although DA are the dominant spectral 
type of white dwarfs. The reason for the delay was observa- 
tional bias. Because of the much higher opacity of hydrogen at 
low temperatures around 6000-20000 K, metal spectral lines are 
far weaker in a hydrogen atmosphere at the same abundances. 
Large telescopes were therefore needed to d e tect DAZs in sig- 
nificant numbers ( Zuckerman & Reidl 119981 : IZuckerman et al.L 
l2003HKoesteretalil2005h . 

Since diffusion timescales are always short compared to evo- 
lution timescales, the metals must be supplied from a source 
from outside, with the possible exception of carbo n in the DQs, 
which can be dredged-up from deeper layers (.Koester et al.L 
|1982|) . Assuming the interstellar matter to be the outer source, 
this so-called accretion/diffusion scenario was discussed in 
great d etail in three fundamental papers by Dupuis et al. ( 199^ 
Il993alfbl) . However, there are some severe problems with this 



scenario: the apparent lack or at least significant underabundance 
of hydrogen in the accreted matter, and the lack of any correla- 
tion between the location of the white dwarfs and the conditions 
of the ISM. Alternative scen arios were therefor e discussed, such 
as the accretion of comets (lAlcock et al.L Il986h . or a tidally dis- 
rupted asteroid (lJurall2003h . A thorough discus s ion o f these al- 
ternative models was given bv lZuckerman et aTl (l2003h . 

These objects have become the subject of renewed interest 
following the detection of infrared excess radiation indicative of 
the presence of circumstellar dust ("Zuckerman & Becklin!, 1987 
Becklin et al., 2005; Kilicetal., 2005, 2006; von Hippel etal 
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120071: lJura et al.L l2007i) . Although the spectra from optical to 
IR can be fairly well modeled with simple disk models (e.g., 
lJura et all l2007h . the exact geometrical distribution was only 
confirmed with the detection of gaseous metal disks around hot- 
ter white dwarfs, which clearly show the signature of Keplerian 
rotation (iGansicke et al.ll20()6f) . 

For all stars exhibiting metal traces in their atmospheres and 
infrared excesses due to circumstellar dust, the source of the cur- 
rent accretion is obviously the matter suiTounding the star. The 
most plausible origin - and the only one discussed in the cur- 
rent literature - is an asteroid remaining from a former planetary 
system, which has been tidally disrupted by the white dwarf. If 
metals are found in the atmosp here , but no infrared excess, the 
situation is unclear However, lJural (|2008|) argues that in such 
a case the circumstellar material may be provided by a higher 
number of small asteroids instead of a single large one. 

One should not forget, however, that the accretion process is 
not well understood. We are not aware of a study of accretion 
from the ISM for realistic conditions in the solar neighborhood, 
with a mixture of gas and dust, and different phases of the ISM, 
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but observationally we know that on all scales in the universe, ac- 
cretion is always accompanied by the formation of an accretion 
disk. We should therefore not exclude from our current consid- 
erations the possibility of some kind of compromise between the 
current two models, accretion from the ISM and accretion from 
the debris disk of a former asteroid. 

It is widely accepted that accretion from the ISM occurs for 
almost all white dwarfs during their lifetimes because more than 
half of all DB (and with higher quality observations perhaps all) 
show traces of hydrogen, whi ch according to o ur current under- 
standing can only be accreted jVoss et al.l 20071) . This is also true 
for the significant amount of hydrogen detected in GD 362. If the 
accreted metals indeed originate in a disrupted asteroid, one has 
to assume an independent source for the hydro gen, unless the as - 
teroid contained a substantial amount of water ( Jura et alll2009l) . 
The approximately solar ratio of metals to hydrogen (see below) 
is in both scenarios purely accidental. 

An important tool in understanding the origin of this dust 
is the determination of its composition by the analysis of the 
atmospheric abundances. The latter is straightforward (if high- 
resolution, high signal-to-noise spectra can be obtained!), but 
the connection between atmospheric and accreted composition 
requires a detailed understanding of the accretion and diffusion 
processes. The composition of the accreting matter need not nec- 
essarily be identical to that of the circumstellar matter (see e.g., 
[Alcock & Illarionov, 1980). And the composition observed in 
the stellar atmosphere is modified by diffusion out of the outer 
reservoir, which occurs at different rates (diffusion timescales) 
for different elements. This is the par t of the problem studied i n 
this paper, which extends the work of Koe ster & Wilkenl (l2006h . 

Since the study of Koester & WiUcen ( 200^ infrared excess 
radiation has been identified in many more hydrogen-rich and 
helium-rich white dwarfs with metals. Improved observations 
have allowed us to identify 15 heavy elements and place up - 
per limits on a few more in GD 362 (Zucker man et al.L l2007l) . 
For many of those elements, diffusion timescales are unavail- 
able in the literature, inhibiting any firm conclusions about the 
composition of the circumstellar matter. We therefore calcu- 
lated atmosphere and envelope models for hydrogen and he- 
lium atmosphere white dwarfs throughout the interesting tem- 
perature range and calculated the timescales for many ele- 
ments. In this work, we follow the ground work laid by the 
Montreal group, using their data on the Coulomb collision in- 
tegrals (^aquette et al., 1986a) and to a large extent the methods 
outlined in lPaquette et all (Il986bl) . 

2. Envelope models 

To determine diffusion timescales, we need to model the struc- 
ture of the outer layers of the white dwarfs, including the 
complete convection zone. The outer boundary is determined 
from atmospheric models, which assume Local Thermodynamic 
Equilibrium (LTE), hydrostatic equilibrium, convective energy 
transport, and a detailed calculation of the radiative energy trans- 
fer. An up-to-date desc ription of i nput p hysics, data, and numeri- 
cal methods is given bv lKoesterl (|2009|) . The parameters of these 
models are effective temperature Tgff and surface gravity log^, 
in addition to the element abundances. The models are calcu- 
lated down into the star to very large Rosseland optical depth 
tr, usually between 1000 and 1500. 

The value of physical parameters such as pressure, tempera- 
ture, and density, at some defined value tr, which is a free pa- 
rameter chosen to be between 1 and 500, is the starting point for 
the integration of the deeper layers. In this integration, we also 



need to know the mass and radius of the star, which w ere ob- 
tained from the finite temperature mass-radius relations of lWoodI 
(1995). 

We use the standard equations of stellar structure for the con- 
servation of mass, hydrostatic equilibrium, and energy transfer. 
The fourth equation, of energy conservation, is replaced by the 
assumption that the energy flux l{r) at the radius r is proportional 
to the mass m{r) inside that radius 



l(r) — — m(r), 
M 



(1) 



for total mass M and luminosity L. This is almost equivalent to 
assuming a constant luminosity throughout the outer layers. The 
equations are integrated inward from the boundary condition un- 
til the fractional mass in the envelope AM/M = (1 - m(r)/M) 
reaches approximately 10 The algorithms and variable defini- 
tions, as well as the program code, are taken in large part from 
our white dwarf structure code, which has been used in numer- 
ous projects. Updated input physics included in our code, are the 
following: 



- the equation of state (EOS) is that of ISaumon et al. I (ll995h . 
which is adapted from the tables obtained from the AAS CD- 
ROM series Vol. 5 (1995). This is probably the most sophis- 
ticated EOS for hydrogen/helium mixtures available today. 

- absorption coefficients we re obtained from the OPAL project 

esias & Rogersl Il996h. supplemented at low tempera- 



tures with the tables of Ferguson et al. (2005^ 



- thermal conductivity data are from Potekhin et al.l (119991) : 
the tables were obtained from their websiteQ 



The OPAL tables were slightly extended at high densities 
by data calculated for H/He mixtures with the help of the Los 
Alamos TOPS program on their web-siteQ Nevertheless, in 
some of our envelope models, the structure falls into a regime 
of the opacity tables, where the assumptions become invalid and 
the data are unreliable. In this case, the highest reliable density 
value at the given temperature was used in the extrapolation. 
Fortunately this had no effect, since in this regime the opacity 
is always high and the models are convective, the temperature 
gradient being practically adiabatic. We confirr ned this findin; 
bousing instead old Los Alamos opacity tables ( iCox & Stewarl 
ll970l) . and no significant changes to the models were found. 

2. 1. The convection zone 



The convection zone in white dwarfs - both the atmosphere and 
envelope - is calculated with a version of the mixing-length ap- 
proximation. It has become customary to denote the specific 
version as, e.g., ML2/q' - 0.6 or in shorthand as ML2/0.6. 
Here ML2 represents a choice of three di mensionless parame- 
ters , which are defined and explained in .Fontaine et al.l (Il98ll) 
and lTassoul et al.l (Il990l) . the number (0.6) is the mixing - length 
as a multiple of the press ure scale height. iKoester et al.l (1 19941) 
and lBergeron et al] (1 19951) demonstrated that a version with "in- 
termediate" efficiency of energy transport is the best choice for 
DAs in ensuring a consistent fit to optical and UV spectra simul- 
taneously. The Bergeron choice ML2/0.6 is the de-facto standard 
for DA model atmospheres today, and our "standard" grid in this 
paper also uses this version in atmosphere and envelope calcula- 
tions. 



http://www.iofFe.rssi.ru/astro/conduct 
|http://www.t4.1anl.gov/cgi-bin/opacity/astro.pI| 
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To our knowledge there is no comp arable published s t udy of 
DB white dwarfs. However, there is in Beauch amp et all d 19991) 
a reference to spectroscopic fits of optical and UV data, which 
favors the version ML2/1.25. We note that these conclusions for 
DAs and DBs are based on observed spectra and thus describe 
the atmospheric layers, which produce the emerging light, i.e., 
above tr 1 . 

It is well established that the deeper layers of the atmosphere 
and envelope are not necessarily correctly described by the same 
version of the mixing-length approximation. The atmospheric 
pai-ameters of G 29-38 are T^tf = 11 700 K, logg = 8.10, when 
analyzed with ML2/0.6 models. The convection zone in such 
a model is extremely thin, and has a lower boundary at around 
optical depth 10. Such a thin cvz has a thermal timescale, defined 



Jq cpTdm 



(2) 



of less than Is. However, G 29-38 is a variable ZZ Ceti star, with 
a shortest period a little larger than 200s. Winget et al. ( 1983) 
and iTassoul et al.l (fT99 0) argued that the thermal timescale of 
the cvz should be similar to the period, that is the cvz should be 
much deeper than indicated by the atmospheric ML parameters. 
The same conclusion is reached by fitting the non-linear light- 
curve ,]Jh£jherm^timesc at 11 700 K is predicted to be a: 
100s (lMontgomervll2005h . 

From a completely different point of v iew, a similar con- 
clusion was reached bv iLudwig et al.l (Il994l) . By comparing the 
mean temperature structure in two-dimensional hydrodynamic 
simulations of the outer layers of a DA white dwarf with model 
atmospheres using MLT, they found that no single MLT version 
can describe the entire structure of the cvz. The efficiency, or in 
the MLT parameterization the mixing-length, must increase with 
depth. 

We therefore calculated a second set of DA envelopes, where 
the MLT version was switched to ML2/2.0 in the envelope cal- 
culation. The result depends strongly on the exact structure of 
the atmosphere model and the location of the "matching point", 
the starting point in the envelope integration. If this layer is be- 
low the atmospheric cvz, e.g., deeper than tr k 10 in the case 
of the G29-38 model, or at such large depth that the convection 
is nearly adiabatic, the envelope integration will not produce any 
deepening of the cvz, regardless of the efficiency of the MLT ver- 
sion. The matching must occur within the super-adiabatic part of 
the atmospheric convection zone for the change in efficiency to 
influence the model structure. We considered the optical depth 
of this layer as a free parameter and calibrated it by demanding 
a thermal timescale of ^ 100s for the G29-38 model. This is of 
course unsatisfactory, and in the future we hope to find a more 
consistent description with a smooth change of efficiency from 
shallow to deeper layers. Nevertheless, our choice is supported 
by the pulsation properties of the DAs and should provide more 
realistic estimates of diffusion timescales. 

T he situation is more favorabl e for the DB stars than for the 
DAs. lBenvenuto & AlthausI (Il997h compared results for thermal 
timescales and cvz depths i n DBs between the m ore so phisti- 
cate d convection theory of ICanuto & Mazzitelhl (ll99lL 1 19921) 
and ICanuto et alj (11996.) . with different simple MLT versions. 
They concluded that a convective efficiency between ML2/1.0 
and ML2/2.0 reproduces the results and predicts the correct lo- 
cation of the blue edge of t he DB instability strip. This was 
confirmed by ICorsico et alj (l2008l) . who concluded that only 
ML2/1 .25 predicts the correct location. A similar result was also 



obtained from the light-curve fitti ng of the prototype variable DB 
star GD 358 ( Mon tgomervl l200"7h . We thus chose this version for 
both the atmosphere and the envelope calculations. 

More fundamentally, it is well known that the MLT ap- 
proximation provides a poor description of the various as- 
pects of a real convection zone. For the case of DA white 
dwarfs, this was studied with an extensive comparison of two- 
dimensi onal radiation-hydro dynamic simulations with ID struc- 
tures bv iFrevtag et al.l (11996) . Even the definition of the lower 
boundary is ambiguous: depending on whether one uses the clas- 
sical stability criterion, the layers with significant convective 
flux, or the layers with non-zero velocities, the resulting mass in 
the "convection zone" can differ by orders of magnitude. While 
the temperature structure (related to convective flux) is proba- 
bly the most important quantity for the pulsational properties, 
for diffusion timescales, the mixed region with non-zero veloc - 
ity is relevant. In the example studied by Frevtag et al.l (Il996l) . 
the mass in the latter is 300 times higher than in the former. 

The general result that the velocity field extends far below 
the lower limit of the unstable region and even the flux over- 
shoot was also confirmed by Montgomery & Kupka (2004) with 
their non-local model of convection in DAs. This implies that 
the diffusion timescales in stars with convection zones may be 
orders of magnitude larger than estimated with our current MLT 
approximations. 



3. Diffusion coefficients and timescales 

For the calculation of diffusion timescales w e fol low closely the 
fundamental works of iPaquette et al.l ( Il986ah and lPaquette et al.l 
('1986b'). From the tables in the first paper, we take the fit 
coefficients for the calculation of Coulomb collision integrals 
and the diffusion coefficients as well as thermal diffusion co- 
efficients. The equations to calculate diffusion velocities are 
taken from the second paper with two minor modifications. 
First, we use what the authors call the "second method" for 
the thermal diffusion coefficien t. This is based on the approach 
by IChapman & CowlingI (Il970l) . and the necessary data can be 
found in the first paper cited above. The second change is purely 
cosmetic. To determine the effective charge of the trace element 
2, we calculate an effective ionization potential XsS as described 
in their Eqs. 19-21, and compare this with the true ionization 
potentials ;t'(Z) of the ions of element 2 with charge Z. Rather 
than taking the effective charge as Z, if XesiT^lxi^) ^ 1 '^"'^ 
Xis{2 + \)lx{'Z, H- 1) > 1, we interpolate a non-integer Zgff be- 
tween Z and Z -H 1 . Diffusion coefficients and diffusion velocities 
are then calculated for ions Z and Z + 1 and a weighted aver- 
age taken depending on the value of Zgff. The only purpose of 
this approach is to avoid unphysical wiggles and steps in the re- 
lation between diffusion timescales and effective temperatures 
(see e.g.. Fig. 4 in Paquette et alj (fl986 b)). 

The diffusion vel ocity Vdis is obtained using Eq. 4 in 
iPaquette et alj ( Il986bl) . neglecting the concentration term, since 
we consider only diffusion of trace elements. The diffusion time 
scale is 



TdifF 



(3) 



where r is the local radius at the bottom of the convection zone 
and p is the local mass density. Results for the stellar models in 
Tables [1] to [3] are given in Tables |4] to |6] for six important ele- 
ments; data for other models and/or other elements can be re- 
quested from the author. 
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Table 1. Conditions at the bottom of the convection zones in DA 
models with surface gravity \ogg = 8 and the standard mixing 
length version ML2/0.6 for atmosphere and envelope. The sec- 
ond column is the fractional mass q = logM^vz/M in the zone 
or down to an optical depth tr = 5, whichever is deeper The 
next columns give the logarithms of gas pressure, temperature, 
and mass density at this level. All physical quantities are in cgs 
units. 





logM„z/M 


logP 


logr 


logp 


6000. 


-7 ^70 


14. SIO 


S 014 


^0^ 


7000. 


-8.423 


13.655 


5.758 


-0.293 


8000 


-Q 007 


1 ^ 076 


S 66S 


-0 77^ 


9000. 


-9.719 


12.359 


5.543 


-1.371 


QSOO 


-10.278 


1 1 7Q0 


S 4^8 


-1 896 


10000. 


-11. 063 


1 1.014 


5.289 


-2.466 


10'SOO 


- 1 7 004. 




S 008 


^10 

J.J i\j 


11000. 


-13.392 


8.684 


4.870 


-4.386 


1 1 100 


1 J.KJO i 


8 ^QS 


4.821 


-4 670 


1 1200. 


-13.975 


8.101 


4.771 


-4.876 


11 300. 


-14.288 


7.788 


4.718 


-5.136 


11400. 


-14.638 


7.438 


4.653 


-5.422 


11500. 


-15.049 


7.026 


4.559 


-5.739 


1 1600 


-15.427 


6 640 




-6 01^ 


1 1700. 


-15.675 


6.401 


4.393 


-6.196 


1 1 xoo 


-1 S 8'SO 

i J.O 


6 775 




-6 ^^0 


11900. 


-15.947 


6.129 


4.332 


-6.404 


12000. 


-16.006 


6.070 


4.321 


-6.450 


12500. 


-16.065 


6.010 


4.316 


-6.504 


13000. 


-16.070 


6.006 


4.326 


-6.521 


14000. 


-16.009 


6.067 


4.361 


-6.500 


15000. 


-15.949 


6.126 


4.391 


-6.472 


16000. 


-15.891 


6.185 


4.420 


-6.444 


17000. 


-15.839 


6.237 


4.447 


-6.420 


18000. 


-15.794 


6.282 


4.473 


-6.400 


19000. 


-15.754 


6.322 


4.496 


-6.385 


20000. 


-15.719 


6.357 


4.518 


-6.371 


21000. 


-15.687 


6.389 


4.538 


-6.361 


22000. 


-15.658 


6.417 


4.558 


-6.351 


23000. 


-15.632 


6.444 


4.575 


-6.343 


24000. 


-15.607 


6.469 


4.592 


-6.335 


25000. 


-15.584 


6.492 


4.609 


-6.328 



4. Elementary considerations for the 
accretion/diffusion scenario 

The possibility of this scenario providing an explanation of 
the observed metals in cool white dwarfs was studied exten- 
sively in a se r ies of three pap ers by Dupuis and collaborators 
jDupuis et all 1 19921 Il993ailbl) . In the first two of those papers, 
the authors solved the complete problem of time-dependent ac- 
cretion and diffusion, by considering the spatial and temporal 
distribution of a heavy element during the white dwarf cooling 
evolution. They also showed that for trace elements in the outer 
convection zone - which is always homogeneously mixed - the 
computation can be replaced by a much simpler local calculation 
at the bottom of the convection zone. We use this approximation 
here to derive some basic consequences of debris disk accretion 
onto white dwarfs, which are not always appreciated in the cur- 
rent literature. 

The basic equation for the mass abundance Xcniet) of 
heavy element el in the convection zone becomes (see Eq. 1 in 



Table 2. Similar to Table[Tl but with a more efficient convective 
energy transport in the deeper layers, as explained in the text. 
For Teff > 13000K the entries are identical with Table[T|and are 
not repeated here. 





logM„z/M 


logP 


logr 


logp 


6000. 


-7.564 


14.516 


5.916 


0.398 


7000. 


-8.376 


13.703 


5.769 


-0.256 


8000. 


-8.876 


13.202 


5.691 


-0.673 


9000. 


-9.429 


12.649 


5.602 


-1.140 


9500. 


-9.803 


12.275 


5.536 


-1.448 


10000. 


-10.328 


11.749 


5.437 


-1.877 


10500. 


-11.040 


11.037 


5.303 


-2.457 


11000. 


-11.945 


10.132 


5.135 


-3.198 


11100. 


-12.151 


9.926 


5.097 


-3.367 


11200. 


-12.364 


9.712 


5.058 


-3.543 


11300. 


-12.587 


9.490 


5.019 


-3.727 


11400. 


-12.819 


9.258 


4.978 


-3.919 


11500. 


-13.058 


9.018 


4.936 


-4.118 


11600. 


-13.301 


8.775 


4.895 


-4.320 


11700. 


-13.544 


8.532 


4.853 


-4.522 


11800. 


-13.806 


8.270 


4.809 


-4.743 


11900. 


-13.996 


8.080 


4.777 


-4.902 


12000. 


-14.229 


7.847 


4.738 


-5.097 


12500. 


-15.555 


6.521 


4.437 


-6.122 


13000. 


-16.070 


6.006 


4.326 


-6.521 



Table 3. Similar to Table[Tl but for DB white dwarfs. Convection 
parameters are ML2/1.25 for atmosphere and envelope models. 





log M„z/M 


logP 


logr 


logp 


6000. 


-4.840 


17.243 


5.938 


2.868 


7000. 


-4.754 


17.330 


6.098 


2.909 


8000. 


-4.760 


17.324 


6.275 


2.882 


9000. 


-4.789 


17.296 


6.440 


2.825 


10000. 


-4.833 


17.253 


6.551 


2.754 


11000. 


-4.954 


17.131 


6.598 


2.636 


12000. 


-5.132 


16.954 


6.603 


2.486 


13000. 


-5.332 


16.752 


6.589 


2.322 


14000. 


-5.545 


16.538 


6.565 


2.151 


15000. 


-5.780 


16.303 


6.532 


1.964 


16000. 


-6.059 


16.023 


6.486 


1.743 


17000. 


-6.413 


15.668 


6.421 


1.462 


18000. 


-6.884 


15.196 


6.331 


1.086 


19000. 


-7.572 


14.507 


6.202 


0.531 


20000. 


-8.657 


13.421 


5.992 


-0.333 


21000. 


-10.046 


12.031 


5.727 


-1.470 


22000. 


-11.368 


10.708 


5.493 


-2.567 


23000. 


-11.581 


10.495 


5.463 


-2.751 


24000. 


-11.746 


10.330 


5.441 


-2.895 


25000. 


-11.911 


10.165 


5.419 


-3.038 



iDupuis et al.L 1 19921 with slightly different notation) 

Mcvz = M(el) - Anr^ X,M t)pv(el) 

at 

. M(.o - (4) 

Tdiff 

where the first term on the right side is the accretion rate of that 
element, and the second term is the rate of gravitational settling 
(where we use a positive velocity to represent elements diffusing 



D. Koester: Diffusion in white dwarfs 



5 



Table 4. Diffusion timescales at the base of the convection zone 
(or at tr = 5) for DA models with log^ = 8, using the standard 
assumption of ML2/0.6 for atmosphere and envelope models. 
Columns 2-6 give the logarithm of the timescales in years for 
the elements C, Na, Mg, Si, Ca, and Fe. 





C 


Na 


Mg 


Si 


Ca 


Fe 


6000 


4.50 


4.34 


4.29 


4.21 


4.14 


4.00 


7000.0 


3.88 


3.60 


3.54 


3.52 


3.45 


3.33 


8000 


3.49 


3.16 


3.07 


3.12 


3.06 


2.95 


9000.0 


3.03 


2.70 


2.60 


2.66 


2.61 


2.50 


9500.0 


2.68 


2.31 


2.21 


2.30 


2.24 


2.13 


10000.0 


2.18 


1.74 


1.65 


1.79 


1.68 


1.59 


10500 


1.55 


0.98 


0.95 


1.15 


0.94 


0.88 


11000.0 


0.55 


-0.00 


0.05 


0.32 


-0.04 


-0.03 


11100.0 


0.33 


-0.22 


-0.13 


0.14 


-0.25 


-0.24 


11200.0 


0.10 


-0.46 


-0.33 


-0.07 


-0.46 


-0.44 


11300 


-0.14 


-0.72 


-0.55 


-0.33 


-0.69 


-0.67 


11400.0 


-0.41 


-1.01 


-0.79 


-0.61 


-0.94 


-0.92 


11500.0 


-0.74 


-1.35 


-1.07 


-0.98 


-1.25 


-1.22 


11600.0 


-1.03 


-1.64 


-1.32 


-1.31 


-1.51 


-1.57 


11700.0 


-1.35 


-1.84 


-1.50 


-1.52 


-1.70 


-1.80 


11800.0 


-1.54 


-1.97 


-1.63 


-1.66 


-1.83 


-1.95 


11900.0 


-1.65 


-2.05 


-1.70 


-1.74 


-1.90 


-2.04 


12000.0 


-1.70 


-2.10 


-1.74 


-1.79 


-1.95 


-2.09 


12500 


-1.74 


-2.14 


-1.78 


-1.83 


-1.99 


-2.12 


13000.0 


-1.73 


-2.14 


-1.79 


-1.83 


-1.99 


-2.12 


14000.0 


-1.61 


-2.08 


-1.73 


-1.76 


-1.93 


-2.04 


15000.0 


-1.48 


-2.03 


-1.69 


-1.70 


-1.89 


-1.97 


16000.0 


-1.36 


-1.98 


-1.65 


-1.64 


-1.84 


-1.90 


17000.0 


-1.32 


-1.94 


-1.62 


-1.58 


-1.80 


-1.83 


18000.0 


-1.28 


-1.90 


-1.59 


-1.53 


-1.77 


-1.77 


19000.0 


-1.24 


-1.86 


-1.57 


-1.48 


-1.74 


-1.71 


20000.0 


-1.21 


-1.83 


-1.55 


-1.43 


-1.72 


-1.68 


21000.0 


-1.18 


-1.80 


-1.53 


-1.38 


-1.69 


-1.66 


22000.0 


-1.15 


-1.77 


-1.51 


-1.34 


-1.67 


-1.64 


23000.0 


-1.12 


-1.74 


-1.50 


-1.31 


-1.65 


-1.61 


24000.0 


-1.09 


-1.70 


-1.48 


-1.27 


-1.63 


-1.59 


25000.0 


-1.06 


-1.67 


-1.47 


-1.24 


-1.60 


-1.57 



downward). Assuming that t, M, v are constant, the solution of 
this equation is 



-tlriel) + T{el)M{el) r 



-iItM) 



\ (5) 



where the first term on the right side is the starting value of the 
element abundance at the beginning of the accretion phase. 

Some conclusions can be drawn immediately from this sim- 
ple equation. If the accretion rate is constant for some multi- 
ple of the diffusion timescale 5 for practical purposes, given 
the typical uncertainties in observed elemental abundances), the 
abundance will approach an asymptotic ("steady state") value 
(we omit the index cvz from X for simplicity) of 



X{el, oo) = 



T{el) M(el) 



and for the ratio of elements 1 and 2, we obtain 

Xiell) _ T{ell)M{ell) _ T(eZl) X^^^iell) 
X{el2) ~ T(el2)M(el2) ~ T{el2) Xacc(eG)' 



(6) 



(7) 



where in the last step, we replaced the accretion rate by the ele- 
ment abundances in the accreted matter If the diffusion times are 
less than a few years, we can reasonably assume that this steady 



Table 5. Similar to Table HI but with more efficient convection 
(ML2/2.0) below tr = 2, leading to deeper convection zones in 
closer agreement with expectations for the variable ZZ Cetis. 





c 


Na 


Mg 


Si 


Ca 


Fe 


6000.0 


4.50 


4.34 


4.29 


4.21 


4.15 


4.01 


7000.0 


3.91 


3.64 


3.58 


3.55 


3.49 


3.37 


8000.0 


3.57 


3.25 


3.16 


3.20 


3.14 


3.03 


9000.0 


3.22 


2.91 


2.81 


2.85 


2.82 


2.70 


9500.0 


2.98 


2.67 


2.56 


2.61 


2.57 


2.46 


10000.0 


2.65 


2.30 


2.20 


2.27 


2.22 


2.11 


10500.0 


2.19 


1.78 


1.68 


1.81 


1.72 


1.62 


11000.0 


1.64 


1.12 


1.05 


1.24 


1.07 


0.99 


11100.0 


1.51 


0.95 


0.91 


1.11 


0.92 


0.85 


11200.0 


1.37 


0.79 


0.77 


0.98 


0.75 


0.70 


11300.0 


1.19 


0.61 


0.61 


0.84 


0.58 


0.54 


1 1400.0 


1.00 


0.43 


0.45 


0.69 


0.39 


0.38 


11500.0 


0.82 


0.25 


0.29 


0.54 


0.22 


0.22 


11600.0 


0.64 


0.07 


0.12 


0.38 


0.03 


0.04 


11700.0 


0.45 


-0.10 


-0.03 


0.23 


-0.14 


-0.13 


11800.0 


0.25 


-0.30 


-0.21 


0.06 


-0.33 


-0.31 


11900.0 


0.10 


-0.46 


-0.34 


-0.07 


-0.47 


-0.45 


12000.0 


-0.08 


-0.66 


-0.50 


-0.26 


-0.64 


-0.62 


12500.0 


-1.15 


-1.74 


-1.42 


-1.41 


-1.61 


-1.68 


13000.0 


-1.73 


-2.14 


-1.79 


-1.83 


-1.99 


-2.12 



Table 6. Similar to TablelH but for non-DA models without outer 
hydrogen layer (DB, DC). 





C 


Na 


Mg 


Si 


Ca 


Fe 


6000.0 


6.50 


6.57 


6.55 


6.59 


6.56 


6.56 


7000.0 


6.48 


6.52 


6.54 


6.52 


6.53 


6.49 


8000.0 


6.40 


6.40 


6.41 


6.40 


6.41 


6.33 


9000.0 


6.34 


6.32 


6.32 


6.32 


6.32 


6.22 


10000.0 


6.31 


6.28 


6.28 


6.29 


6.27 


6.15 


11000.0 


6.25 


6.21 


6.22 


6.20 


6.18 


6.06 


12000.0 


6.17 


6.13 


6.13 


6.09 


6.05 


5.94 


13000.0 


6.06 


6.01 


5.98 


5.96 


5.88 


5.80 


14000.0 


5.96 


5.88 


5.85 


5.85 


5.72 


5.64 


15000.0 


5.84 


5.72 


5.71 


5.71 


5.53 


5.46 


16000.0 


5.68 


5.54 


5.54 


5.49 


5.33 


5.24 


17000.0 


5.51 


5.35 


5.32 


5.23 


5.11 


4.98 


18000.0 


5.26 


5.05 


4.99 


4.89 


4.83 


4.64 


19000.0 


4.73 


4.62 


4.56 


4.44 


4.43 


4.22 


20000.0 


4.06 


3.99 


3.91 


3.77 


3.78 


3.56 


21000.0 


3.22 


3.12 


3.00 


2.86 


2.96 


2.68 


22000.0 


2.36 


2.15 


2.05 


1.89 


2.03 


1.82 


23000.0 


2.21 


2.00 


1.90 


1.74 


1.89 


1.67 


24000.0 


2.07 


1.85 


1.75 


1.60 


1.75 


1.53 


25000.0 


1.90 


1.66 


1.57 


1.42 


1.57 


1.37 



state has been reached. Since we can calculate the timescales 
from the parameters of the star, we can thus determine the abun- 
dances in the accreted matter from the observed abundances in 
the star. Even more simply, since the timescales for the com- 
monly observed elements of Ca, Mg, and Fe are often within a 
factor of two, we can assume the observed abundances (which 
are rarely more accurate than a factor of 2) to be a first approxi- 
mation to the accreted abundances. 

However, all of these conclusions depend critically on the as- 
sumption that the steady state is reached, and are therefore valid 
only in the case of short diffusion timescales. It is instructive to 
consider by a simple example, the consequences of being unable 
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Fig. 1. Simple example for the accretion/diffusion scenario with 
two elements. 



to assume steady state (i.e., in all cases where the timescales ex- 
ceed a few years). 

For this exercise, we studied the abundances of two elements, 
where the diffusion timescales differed by a factor of two. We 
started with zero abundance and switched accretion on for 10 
times the shorter timescale. Then accretion is switched off again. 

We were able to distinguish three phases (see Fig.[T]| 

1. for t <K Ti, T2, expansion of the exponentials in the equation 
above shows that the diffusion timescales cancel; the element 
abundances increase linearly with time, with the slope given 
by the accretion rates (or the abundances in the accreted mat- 
ter) 

2. for t X T, the abundances approach the steady state values, 
which are reached after w 5t. In this phase, the abundances 
of the accreted material are modified by the ratio of the dif- 
fusion timescales, which is here equal to 2 

3. when the accretion is switched off, the abundances de- 
crease exponentially. Because of the difference in diffusion 
timescales, the ratio also increases exponentially and can 
reach very large values - in this simple case > 100 after 10 
diffusion timescales. At this point, the element abundances 
are only 2-4 orders of magnitude below their maximum val- 
ues, which would be observable in many cases. 

The third phase would also describe a different scenario, 
where accretion occurs at a high rate in a short time, such as the 
accretion of an entire asteroid at once. The abundances would 
reflect exactly the accreted abundances during this time and start 
the exponential decay from there. In such a case we would be 
unable able to determine the composition of the accreted matter. 

We note that the above example is not extreme. In particular, 
when including light or very heavy elements in the comparison 
the factor between timescales can be a 4. And we also do not 
need to consider the extreme case of totally switching accretion 
on and off. Even a change in the accretion rate, for example by 
a factor of 100 leads to an intermediate change in the abundance 
ratio of a factor > 10, before the asymptotic ratio of 2 is reached 
again. 

The conclusion of this simple exercise is that it is only possi- 
ble to infer abundances of the accreted matter if we can safely as- 
sume that the accretion rate has been constant for several times 
the diffusion timescales of all elements involved, and that there- 
fore the steady state has been reached. For diffusion timescales 
longer than a few decades at most this can never be assumed. 



Equating the observed atmospheric abundance ratios to those of 
the accreted matter can in such cases be incorrect by orders of 
magnitude. 

5. Specific examples 

We apply the general results to two specific objects with circum- 
stellar material and atmospheric metal traces, which have been 
discussed in numerous recent studies: G 29-38 and GD 362. 



5.1. G 29-38 



Based on three independent deter minations jLiebert et al.L 120051 : 
iKoester et all [19971 l2006h . we use the atmospheric pa- 
rameters Teff = 11 690 K, log ^ = 8.11. An atmosphere/envelope 
model was calculated with these data and Table[7]c ollects the dif- 
fusion times cales and observe d abundances from iKoester et al.l 
(1 19971) and IZuckerman et al.l (|2003|) . Unfortunately, there is 
some confusion about the Ca abundance in the literature - even 
using only studies with contributions from this author The rea- 
son is that the Vienna Atomic Line Database (VALD) contains 
Stark broadening data for the Can H and K resonance lines, 
which differ by about a factor of 10 from most of the Kurucz 
lists. Depending on which data are used, the abundance differs 
by about 0.5 dex. Using simple estimates for Stark broadening 
indicates that the larger broadening constants are more likely 
correct and thus the lower Ca abundance is preferred. 

For the model, we used the version with the more efficient 
convection, which provides a reasonable thermal timescale of 
about 100 s. Since the diffusion timescales are still smaller than 
one year, the conclusions are not qualitatively changed compared 
to the standard version, which was used in IKoester & Wilkenl 
("2006^ There is a significant caveat, however: as discussed in a 
previous section, even though the thermal timescale agrees with 
the pulsation properties, the possible overshooting at the bottom 
might increase the diffusion timescales significantly. Another 
effect in the same direction is the pulsation itself. During the 
pulsational decrease in the effective temperature (up to 500 K; 
Montgomery, 2005), the depth of the convection zone increases 
significantly, leading to a larger mixed zone and an increase in 
the diffusion timescale by almost a factor of 10 compared to the 
equilibrium model. 

Metals were first detected in G 29-38 in 1997 (K oester etall 

Il997i) . Von Hippel & Thompson (2007) claimed changes 
had occurred in the equivalent width of the Can reso- 
nan ce lines during the p a st dec ade, which was questioned by 
Deb es & Lopez-MoralesI (|2008|) . We have seven high-quality 
high-resolution spectra from the ESO VLT and the Keck tele- 
scope, w hich show no significant variati on between 1997 and 
2000; the I von Hippel & ThompsonI (|2007|) result may have been 
influenced by the inhomogeneity of their data, which included 
time-resolved spectroscopy of very low S/N. 

In view of the short diffusion timescales, it seems thus rea- 
sonable to assume that this object at the moment is experienc- 
ing steady state accretion. Because of the similar timescales, the 
predicted accretion abundances do not differ too much from the 
observed values. Taking into account that there is a factor of 3 
diffe rence in the Mg abundance between iKoester et al.l (Il997l) 
and Zuckerman et al.l (|2003|) and a similar uncertainty in the Fe 
abundance (high vs. low dispersion spectra, optical vs. UV spec- 
tra), we can conclude that within the rather large uncertainties, 
the abundances of Ca, Fe, and Mg in the accreted matter are 
compatible with solar abun dances, take n from the compilation 
in Astrophysical Quantities (Coxl |2000|) . 
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Table 7. Diffusion data for the variable ZZ Ceti star G 29-38. 
Column 2 gives the time scale in years, column 3 the observed 
number abundances relative to Ca, column 4 the observed abun- 
dances by mass, also relative to Ca. Column 5 is the predicted 
abundance in the accreting material, assuming steady state be- 
tween accretion and diffusion. The final column shows the solar 
abundances on the same scale. All abundances are in decadal 
logarithms. 



el 


r [yrs] 


log ealm 


log X;„m 




logXe 


Ca 


0.79 


0.00 


0.00 


0.00 


0.00 


Mg 


0.97 


1.10 


0.88 


0.79 


1.00 


Fe 


0.80 


0.61 


0.82 


0.81 


1.33 



5.2. GD362 

GD 362 was dis c overed as a massive DA wit h metals (DAZ) by 
iGianninas et al.1 (|2004 . An infrared excess (iKihc et al.L 120051) 
placed it into the new class of DA white dwarfs with debris 
dis ks. However, an a nalysis of high S/N, high-resolution spec- 
tra (IZuckerman et al.L [2007. Z07 henceforth) showed it to be a 
normal mass, but heUum-rich object. We used the stellar pa- 
rameters for an atmosphere/envelope model, and inferred the 
diffusion parameters in Table |8] We found that the fractional 
mass in the convection zone is <7 = -6.58, and, using the mass 
0.73 Mq from Z07, that the total hydrogen mass in the envelope 
is 3.51 X lO"** Mo, and the total mass of Si is 3.90 x 10^' g. 
The combined mass of the clearly identified elements Mg, Al, 
Si, Ca, Fe, and Ni is ^ 1.8 x 10^^ g, in excellent agreement with 
the estimates in Z07. 

Table 8. Diffusion data for GD 362. The meaning of the columns 
is the same as in Table |7] except that all element abundances are 
presented relative to Si. 



el 


T [yrs] 


log fatm 


logXaim 


logX,,, 


logXo 


C 


1.28 10^ 


<0.20 


<-0.17 


<-0.38 


0.64 


N 


1.2210^ 


<1.70 


<1.39 


< 1.20 


0.20 


O 


1.13 10^ 


<0.70 


<0.45 


<0.30 


1.14 


Na 


9.50 10^ 


-1.95 


-2.04 


-2.11 


-1.33 


Mg 


9.44 lO"* 


-0.14 


-0.21 


-0.28 


-0.04 


Al 


8.57 lO"* 


-0.56 


-0.60 


-0.63 


-1.10 


Si 


7.95 10^ 


0.00 


0.00 


0.00 


0.00 


Ca 


6.07 lO"* 


-0.40 


-0.25 


-0.13 


-1.04 


Sc 


5.83 lO"* 


-4.35 


-4.15 


-4.02 


-4.24 


Ti 


5.45 lO'' 


-2.11 


-1.88 


-1.72 


-2.33 


V 


5.20 lO"* 


-2.90 


-2.64 


-2.46 


-3.29 


Cr 


5.05 lO"* 


-1.57 


-1.31 


-1.11 


-1.61 


Mn 


4.77 10^ 


-1.63 


-1.34 


-1.12 


-1.87 


Fe 


4.67 lO"* 


0.19 


0.49 


0.72 


0.29 


Co 


4.43 lO"* 


-2.66 


-2.34 


-2.09 


-2.31 


Ni 


4.45 lO"* 


-1.23 


-0.91 


-0.66 


-0.98 


Cu 


4.25 lO"* 


-3.34 


-2.99 


-2.72 


-2.98 


Sr 


3.14 lO"* 


-4.58 


-4.09 


-3.69 


-4.15 



Given that the individual abundances vary over more than 
four orders of magnitude, it is remarkable, how close the inferred 
accretion abundances (column 5) are to the solar abundance ra- 
tios. The error for Si was given as 0.3 dex in Z07, and for the 
other elements, it ranges from 0. 1 to 0.4 dex. The only elements 



-5.00 






- 










-5.50 










/' 






>< -6.00 
o 


/ Mg 
Si 






-6.50 


Ca 






-7.00 


Fe 
Sr 

1 


1 


1 \ i\ 



2.00E-H05 4.00E-1-05 6.00E-1-05 8.00E-1-05 

time [yrs} 

Fig. 2. Possible accretion/diffusion history for GD362. The el- 
ement symbols are ordered in the same sequence as the corre- 
sponding curve. 

with abundances that differ significantly from the solar value are 
C and O, for which only upper limits are known, and Ca and Na. 
Since 12 of 14 elements have a ratio of their abundance to that of 
Si higher than for solar values might imply that the Si abundance 
in Z07 is slightly too low, producing a systematic shift. 

The possible "abundance history" of GD 362 is shown in 
Fig.|2] which is in a similar spirit as Fig.[r| We use the elements 
C, Mg, Si, Ca, Fe, and Sr as representatives for different diffu- 
sion timescales and abundances. Accretion is switched on at time 
f = 0, when no heavy elements are present in the atmosphere; the 
rate - the same for all elements - is calculated to ensure a steady 
state abundance of silicon of logX(Si) - -5, the observed value. 
The silicon rate necessary for this is 2.4 x 10'^ g per year. A dif- 
ferent rate would shift all curves vertically by the same amount. 
In reality, the elements would of course be accreted with differ- 
ent abundances. In our choice of presentation, the difference be- 
tween the curves shows directly the relative change between the 
accreted and atmospheric abundance ratios. When all the curves 
are approximately in agreement, the observed abundance ratios 
are identical to the accreted ratios. In this specific case, this is 
true for times f ^ 3 x 10"^ yrs. After that time, the steady state is 
approached, and the ratio differs from the accreted ones by the 
ratios of the diffusion timescales, which are factors of 2-4 in this 
case. 

After the accretion is switched off, the picture changes sig- 
nificantly. About 4x10^ yrs after the end of accretion, the ratios 
can differ from the accreted abundance ratios by factors of up 
to 200. This decline is fast on astronomical timescales, but slow 
enough to be invisible over a few decades; we thus do not have 
any direct indication of the phase we are currently observing. It 
should also be noted - as in the simple example above - that a 
change in the accretion rate of a factor of 10-100 would force the 
star into a phase of exponential decline, until a new steady state 
with lower abundances was reached. 

In their s tudy o f accretion from the interstellar matter, 
iDupuis et al.l ( |l993a|) assumed schematically a duration of the 
accretion phase within a dense ISM cloud of lO*" yrs. This would 
correspond to the time required to complete all of the process in- 
dicated in Fig|2] and all phases would be possible for the current 
situation. On the other hand, in the scenario of accretion from the 
debris disk formed by an asteroid (see e.g., Z07 for an extensive 
discussion), the total accreted matter of the identified elements 
alone would be « 7.5 x lO^-' g after 3x10'' yrs, close to the mass 
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l.OOE+05 2.00E+05 3.00E+05 4.00E+05 5.00E+05 6.00E+05 
time [yrs} 

Fig. 3. Alternative history for GD 362, assuming an exponen- 
tially declining accretion rate. Otherwise the meaning of the 
curves is the same as in Fig.|2l 



In this example, the accretion timescale is longer than all 
diffusion timescales. As time increases, the later decline is gov- 
erned by the accretion timescale, which is the same for all ele- 
ments. The abundance ratios therefore do not show the large de- 
viations from the accreted values as in the case, where accretion 
is switched off. Instead of the ratios of the diffusion timescales 
that govern the relative changes in the "steady state", we obtain 
here the ratio 

X{el\,^) ^ T{el\)X.,,,(el\) t,,, - T(el2) 
X(el2, oo) T(el2) X.,,Jel2) r.^c - T(e/1) ' 

which for very large Tacc relative to T{el) recovers the solution 
for the steady state. 

If the accretion timescale is shorter than the diffusion 
timescale of some elements, then at later times for these ele- 
ments the conditions approach an exponential decline after the 
switching off of accretion, and there is the possibility of large 
deviations from the original abundance. This might in the future 
be helpful in providing constraints on the lifetimes of the disks. 



of Ceres, the most massive asteroid in our system. This agrees 
with estimates of the lifetime of the debris disks,_which are of 
the order of 1.5 x 10^ yrs JjuraL [20081: iKilic et al.L 120081) . 

It is thus highly likely that we are witnessing the early phase 
of accretion, and that it will be possible to identify the com- 
position of the accreted matter This view is further supported 
by the visibility of rare elements such as strontium, which have 
short diffusion timescales. Within 10^ yrs of the end of the ac- 
cretion phase, the Sr abundance would decrease by an order of 
magnitude, becoming undetectable. 

We also note that the aforementioned deviations from solar 
abundances for C, N, Na, and Ca cannot be explained by assum- 
ing a post-accretion phase. Deviations would be in the opposite 
direction of expectations: C, O, and Na are underabundant, but 
should increase in this phase, and vice versa for Ca. 

An element of particular importance for the accre- 
tion/diffusion scenario is hydrogen, because as the lightest el- 
ement, it will always remain in the convection zone and not dif- 
fuse downward. Assuming that hydrogen would be accreted with 
solar ratio with respect to silicon, the current hydrogen mass 
present - 6.98 x lO^'^g - would be reached in 16 580 yrs. We 
either are observing GD 362 in the very early phase of accretion, 
or t he accrete d matter is significantly hydrogen-depleted. 

lJura et all (1200 9) considered an exponentially decreasing ac- 
cretion rate, which would be appropriate if the dust disk is not 
replenished from a reservoir farther out in the system. In such 
a case, the diffusion equation can also be solved analytically to 
have the solution (in our notation) 



Xcyziel, t) 



X,MO)e-"''-"'> 

T(el) Moid) Tacc 



r{el) 



[e-'/-- - e-'/^<^')] (8) 



with a timescale for the accretion rate of t^cc and an initial rate 
Mo(eO- We calculated the "abundance histor y" for GD 362 wit h 
the smaller of the two times considered in 'J ura et al ] (I2009h . 
Tacc = 1.5 X 10^ yrs, and show the results in Fig. [3] 

In this case, the abundances never reach an asymptotic state, 
but instead begin to increase, reach a maximum, and then de- 
cline. The maxima are reached at different times for the ele- 
ments, depending on their diffusion timescales. The initial ac- 
cretion rate for Si (and all other elements) had to be increased 
by a factor of two from the case of constant rate, because other- 
wise the Si abundance would never reach the observed value. 



6. Conclusions 

We have extended the calculation of diffusion timescales to 
many more elements than currently available in the literature, to 
provide the data necessary to interpret the observed abundances 
in the increasing number of white dwarfs with traces of met- 
als in their atmospheres and an infrared excess, indicative of the 
presence of a dust disk. With a simple example and the appli- 
cation to two of the most extensively observed objects (G29- 
38 and GD 362), we have demonstrated that the observed abun- 
dances are in general not identical to the accreted abundances. 
If, however, we have good arguments that the steady state be- 
tween accretion and diffusion is reached, the differences are typ- 
ically only a factor of 2-4, in many cases within the errors of 
the determination. This is still true, if the accretion rate changes 
only slowly, on timescales longer than all diffusion timescales 
involved. With knowledge of the timescales, we can improve the 
comparison by calculating the accreted abundances from the ob- 
served ones. With further observations of the highest quality, it 
should thus be possible to resolve the question about the origin 
of the accreted matter. 
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